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ABSTRACT 

We present a model for the formation of large organic molecules in dark clouds. The 
molecules are produced in the high density gas-phase that exists immediately after ice 
mantles are explosively sublimated. The explosions are initiated by the catastrophic 
recombination of trapped atomic hydrogen. 

We propose that, in molecular clouds, the processes of freeze-out onto ice mantles, 
accumulation of radicals, explosion and then rapid (three-body) gas-phase chemistry 
occurs in a cyclic fashion. This can lead to a cumulative molecular enrichment of the 
interstellar medium. 

A model of the time-dependent chemistries, based on this hypothesis, shows that 
significant abundances of large molecular species can be formed, although the complex- 
ity of the species is limited by the short expansion timescale in the gas, immediately 
following mantle explosion. 

We find that this mechanism may be an important source of smaller organic 
species, such as methanol and formaldehyde, as well as precursors to bio-molecule for- 
mation. Most significantly, we predict the gas-phase presence of these larger molecular 
species in quiescent molecular clouds and not just dynamically active regions, such as 
hot cores. As such the mechanism that we propose complements alternative methods 
of large molecule formation, such as those that invoke solid-state chemistry within 
activated ice mantles. 
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1 INTRODUCTION 

A variety of laboratory evidence shows that the catastrophic 
recombination of hydrogen atoms and other accumulated 
radicals in a soli d may abruptly raise the t emperature of the 
solid to ~ 10 3 K. lDulev fc Williams! (|201lD proposed that for 
HAC dust grains in the interstellar medium, such tempera- 
tures are sufficiently high to permit the grains to radiate in 
the so-called Unidentified Infrared Bands (UIBs) associated 
with carb onaceous materials, lying between 3.3 a nd 11 .3 fim. 
Recently, ICecchi-Pestellini. Dulev fc Williamsl (|2012l ) have 
shown that H2 molecules released from HAC grains in such 
abrupt temperature excursions caused by the recombina- 
tion of a sufficient number of accumulated H-atoms could 
not only account for the populations in rotational levels 
of H2 observed in the "hot component" of H2 in diffuse 
clouds, but also produce amounts of CH + and OH compa- 
rable with those observed in diffuse clouds. This theory of 
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radical recombination-dri ven mantle explo sions is very sim- 
ilar to that described by iGreenberd (|l976l ) except that the 
mechanism for mantle explosion is the spontaneous internal 
recombination of trapped hydrogen atoms, rather than an 
external heating source. 

In this paper, we consider an application of the idea 
of abrupt temperature excursions caused by radical recom- 
bination to dark clouds where ices have been deposited 
on dust grain surfaces. If the abrupt t emperature excur- 
sions are of the magnitude discussed bv lDulev fc Williamsl 
i|201ll) then any ices present will be abruptly converted 
to gas and evaporated exp losively. As pointed out by 
ICecchi-Pestellini et alj |2010|), hereafter Paper I, the den- 
sity in the expanding gas can be so high that three-body 
reactions have the potential to create new and more com- 
plex sp ecies from the ice constituents . Here, we extend the 
ideas of lCecchi-Pestellini et all l|2010t ) to the gas formed by 
the exploding ices. 

Interstellar ices are deposited on dust grains in dark 
interstellar clouds in regions where the visual extinction ex- 
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ceeds some critical value, typically a few visual magnitudes. 
The ices in quiescent lines of sight towards low mass stars 
are observed to contain H 2 0, CO, C0 2 , CH 4 , NH 3 , H 2 CO, 
CH3OH, OCN~, a nd some other species of low abundance 
l|6berg et al.ll20l"lh . The ices may also contain radicals cre- 
ated from these molecules by photodissociation, where the 
disso ciating radiation is the c osmic ray-induced radiation 
field l|Prasad fc Tarafdarlll983l l. 

Thus, H2O may give rise to OH, CH4 to CH3, etc., and 
the population of these radicals is associated with the bulk 
of the ice. The various reactions resulting from the accre- 
tion of atoms and molecules from the ambient gas as well 
as processing by ultraviolet photons and cosmic rays must 
always be incomplete and will, for example, lead to the pres- 
ence of trapped H-at oms at weak binding sites on or near the 
surface. According to lDulev fc Williams! (|201ll ). when a crit- 
ical number of these weakly bound H-atoms is accumulated, 
then a localised recombination of a few hydrogen atoms can 
trigger a chemical runaway which releases all of the chem- 
ical energy stored in the grain. This can include both the 
hydrogen recombination energy and the energy stored in 
other radicals as well as the grain substrate. This leads to a 
runaway explosion which abruptly heats the ice mantle and 
grain core to temperatur es on the order of one thousand K. 
iDulev fc Williams! (|201lh show that the number of H-atoms 
required to cause this explosion is equivalent to about 5% 
of the total number of atoms in the grain plus mantle. We 
postulate that the recombination of trapped H atoms to H2 
drives the explosion but that the heavy, and less mobile, 
molecular radicals are released into the gas phase. 

In the work presented here, we follow the time- 
dependent gas-phase chemistry in a dark interstellar cloud; 
we compute the deposition of ices and the accumulation of 
weakly bound H-atoms that are assumed to trigger the ex- 
plosions when their number on a grain reaches the critical 
value. Thus, in this model, the frequency of the episodic 
temperature excursions is controlled by the cloud chem- 
istry. When a temperature excursion occurs, the ices are 
assumed to be instantaneously sublimated, and three-body 
reactions involving radicals produce products in the very 
high density expanding gas. This gas is assumed to undergo 
a free expansion which occurs on a timescale on the order 
of nanoseconds; however, the density in this sublimate is so 
high that many collisions occur before the gas has relaxed 
to more normal interstellar conditions. In this picture, the 
episodic explosions therefore enrich the interstellar gas with 
the products of the three-body reactions, and the material 
undergoing this enrichment is accumulated during the inter- 
val between explosions. The interstellar gas may be enriched 
by a number of successive events, if the timescales for accu- 
mulation of ices and H-atoms (typically, on the order of one 
million years in canonical dark cloud conditions) are suit- 
able. Evidently, if the interval is very short (as will occur 
if the H-atom abundance in the gas is high), little material 
can be accumulated as ices. 

Our model for molecular cloud chem- 
istry has some similarity to that described by 
iGarrod. Widicus Weaver fc Herbstl (2006f) to account 
for chemical complexity in hot cores. The accumulation 
of radicals within ices is adopted in both our model of 
chemistry in molecular clouds and that of Garrod et al. in 
hot cores. In our model, the radicals are released through 



the abrupt warming in episodic explosions, while in the 
Garrod et al. model of hot core chemistry the radicals are 
released in the slow warm-up of material in the early hot 
core phase. 

We describe the model in Section 2. The results are 
presented in Section 3, and we give a discussion of the im- 
plications and make our conclusions in Section 4. 



2 THE MODEL 

To test the hypothesis described above, we first of all con- 
structed a simple 'proof of concept' model of the time- 
dependent chemistry. The first stage considers the chemistry 
(gas-phase, plus freeze-out), for moderately dense conditions 
(density n = 10 3 cm -3 , temperature T = 10 K, extinction 
Av = 3) where we expect the atomic hydrogen abundance 
to be high, so as promote the rapid accumulation of radicals 
in the ice mantles. During this phase limited surface chem- 
istry is allowed to occur, so that full hydrogenation of C, 
N, O, and S-hydrides to CH 4 , NH 3 , H 2 and H 2 S are al- 
lowed. A fraction (10%) of the CO is allowed to react with O 
and OH on the surface of grains to form C0 2 . It is further 
assumed that no (continuous) desorption mechanisms are 
operating. After a specified period (~ 10 4 — 10 6 years) the 
ice mantles (composed of CO, CH4, NH3, N 2 , 2 , H 2 0, Na, 
H 2 CO, C0 2 , HCN, HNC, HNO, H 2 S, C 2 S, HCS, O and OH) 
are released back into the gas phase. No change of tempera- 
ture or density is included and no high-density chemistry is 
included. The cycle is then allowed to repeat. 

The results from this model show that, for most species, 
a limit cycle is achieved after only one or two mantle explo- 
sions. Ions and unsaturated species tend not to be strongly 
affected by the process, but the simple saturated species 
(such as H 2 0, CH4 etc.) are very strongly enhanced. A few 
species, such as CO and C0 2 show a slower, more steady 
rise in abundance over several cycles. 

We have therefore constructed a more realistic and 
comprehensive two-phase model which considers the time- 
dependent chemistry at a single point in a molecular 
cloud. The parameters for this model are given in Ta- 
bic fT] This model utilizes th e LSODE integration package 
l|Hindmarsh fc Petzoldl I1995T ) and has been sub-divided to 
describe the two chemical phases. Phase I represents the 
(standard) dark cloud chemistry, with freeze-out and (lim- 
ited) surface chemistry, whilst Phase II considers the chem- 
istry in the high density, rapidly expanding gas, in the im- 
mediate vicinity of a dust grain following ice mantle subli- 
mation. 

The two phases are physically and chemically distinct 
from each other, but the output from each phase feeds into 
the other as material cycles between the two. The charac- 
teristics of the two phases are described below:- 

2.1 Phase I 

The physical and chemical parameters for Phase I (as given 
in Table QJ are discussed below. 

For the gas-grain interactions, all species are assumed to 
have a sticking coefficient of 1.0 and as with the simple model 
- and for the purpose of clarity - we suppress all continuous 
desorption processes. The freeze-out rates (and hence the H 
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Table 1. Physical parameters in the standard model - see text 
for description 



Parameter Value 



He/H 0.1 

C/H 2.6 x 1CT 4 

N/H 6.1 x 1CT 5 

O/H 4.6 x 10~ 4 

S/H 1.0 x 10~ 7 

Na/H 1.0 x 10~ 7 



Density (rij) 10 4 cm 3 

Temperature (Tj) 10 K 

Ay 3 magnitudes 

Cosmic ray ionization rate (£o) 1.3 X 10~ 17 s — 1 

Initial abundance of H-atoms (rijj q) 1 cm~ 3 (see text) 



H-atom non-recombination probability (ph) 0.1 

Explosion threshold abundance of H (/#) 0.05 

No. of (refractory) atoms per grain (N g ) 10 s 

Mantle radical formation rate (R ra d) 0.01 Myr - 



Average grain radius (a) 0.0083 fim 

Dust surface area per H-nucleon (oh) 8.0 X 10 — 21 cm 2 
Grain albedo 0.5 
CO — > CO2 conversion efficiency (/co 2 ) ".1 



Phase II: Initial density (njj) 10 20 cm -3 

Phase II: Initial Temperature 1000 K 

Phase II: Three-body rate coefficients (/C3s) 10 -28 cm 6 s — 1 



Number of cycles (n cyc .) 5 



created in the ice mantles due to the action of impinging 
cosmic rays. We note that there may be contributions both 
from direct cosmic ray imp act and photolysis by the cosmic 
ray induced radiation field l|Prasad fc Tarafdarll 19831 ). This 
is a cumulative effect and we assume that it is limited to the 
stripping of a single hydrogen atom from saturated species. 
Thus: 

CH4 -> CH 3 + H 
NH3 ->■ NH 2 + H 

etc. This applies to; CH 4 , NH 3 , H 2 0, OH, H 2 CO, H 2 S, 
CH3OH, HCN, HNC, HNO and HCS. The relative abun- 
dance of the radical to the parent saturated species will 
therefore be proportional to the cosmic ray ionization rate 
and the period of exposure, t. This fraction is then given 
by:- 

n-rad R f t \ / Co A 

nsat \ 1 Myr" 1 J V 1.3 x lO" 17 ^ 1 ) 

Note that we have not included the (equally likely) possi- 
bility of the formation of molecular ions in the ice mantles. 
This is a potentially serious omission, but since we have lit- 
tle idea as to the efficiencies and products we have opted not 
to include these processes in our model. However, we also 
note that molecular ions are highly reactive, so the results 
from this study must be regarded as a lower limit to the 
formation efficiencies of complex organic molecules. 

It is further assumed that the hydrogen atoms so- 
generated are released back into the gas-phase and are not 
retained in the ice mantles. These radicals will be crucial 
in determining the efficiency of large molecule formation in 
Phase II, as we postulate that radical-radical three-body re- 
actions will be very much faster (by a factor of > 100 x) 
than reactions between saturated species. We use a value of 
Rrad = 1% per Myear in our standard model. 

In the case of atomic hydrogen we assume that a frac- 
tion (ph) that impacts grains simply sticks as free hydro- 
gen atoms in the ice. The remainder recombines to H 2 
and is returned to the gas phase, as per standard mod- 
els of interstellar c louds. In our standard model we follow 
iDulev fc Williams! (|201lt ) and adopt a value of 0.1 for pu- 
The timescale between mantle explosions is determined by 
the rate of accretion of atomic hydrogen, so the initial H 
abundance is an important parameter. Assuming that chem- 
ical equilibrium initially pertains in the cloud then, for typ- 
ical gas-to-dus t ratios, the abundance of atomic hydrogen is 
given by (e.g. IDulev fc Williamslll984 l 

n H = l( C ° - - ) cm" 3 

V 1.3 x lO-^s- 1 / 

There are 81 gas-phase and 25 solid-state species in the 
dark cloud (Phase I) chemistry which are listed in Table [2] 
Those given in italics are solid-state (frozen-out) species, 
whilst those in bold are the solid-state radicals formed in 
situ. 

The time-dependence of the chemistry is followed un- 
til the atomic hydrogen abundance in the ices reaches some 
threshold value (Jh)- At this point we assume that the hy- 
drogen explosively recombines (with 100% efficiency) to H 2 
and all components of the ice mantles are instantaneously 
heated and fully sublimated. The energy liberated has to 



to H 2 conversion rate) are calculated self-consistently, using 
the mean grain radius (a) and surface area per hydrogen 
nucleon (an) as specified in Table [1] These values are con- 
sistent with standard values for the dust-to-gas ratio an d 
freeze-out rates in dark clouds (e.g. iRawlings et al.lll992h . 
We use a typical value for the cosmic ray ionization rate and 
the reaction network is drawn from the UMIST06 database 
l|Woodall et al.ll2007l ). We also include photodissociations by 
the secondary radiation field gen erated by the cosmic ray 
ionization and excitation of H 2 (|Prasad fc Tarafdarl 119831 ) 
which have an inverse dependence on an and (1— uj), where 
u) is the average grain albedo. 

We assume that full hydration of atoms and simple hy- 
drides to CH4, NH3, H 2 Q and H 2 S occurs on gra in surfaces. 
As in other studies (e.g. IRawlings & Kcto 20121) we assume 
that a fraction (fco 2 ) of the CO that impacts a grain and 
interacts with surface O or OH is converted into C0 2 . Some 
models of hot cores (e.g. IViti fc Williams! 1 1999D invoke the 
partial conversion of CO to CH3OH and/or H 2 CO on the 
surface of grains to explain the high gas-phase abundances 
that are observed. Whether or not such processes are effi- 
cient dWatanabe fc Kouchil2002T) we have not included them 
in our model as we wish to distinguish the role of the pro- 
posed mechanism in the generation of larger organic species. 

The extinction (Av) is taken to be 3 magnitudes for 
ni — 10 4 cm -3 and 10 for higher densities. It is assumed 
that the chemistry is 'dark' - i.e. well outside any photon- 
dominated regions (PDRs), so that the photodissociation 
rates for H 2 and CO and the photoionization rate for C are 
all set to zero. 

During this phase we assume that reactive radicals are 
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Table 2. Chemical species in Phase I 



Table 3. Chemical species in Phase II 



H, H 2 , H+, H~, H+, H+, He, He+, Na, Na+, 
C, C+, C", CO, CO+, CH, CH+, CH 2 , CH+, 
CH 3 , CH+, CH 4 , CH+, CH+, 
N, N+, NH, NH+, NH 2 , NH+, NH 3 , NH+, NH+, N 2 , N+, N 2 H+, 
O, 0+, 2 , 0+, OH, OH+, H 2 0, H 2 0+, H3O+, 
HCO, HCO+, H 2 CO, H 2 CO+, C0 2 , CO+, HCO+, 
CN, CN+, HCN, HCN+, HNC, 
NO, NO+, HNO, HNO+, HCNH+, H 2 NC+, HNCO+, H 2 NO+, 
S, S+, HS, HS+, H 2 S, H 2 S+, H 3 S+, 
CS, CS+, C 2 S, C 2 S+, HC 2 S+, HCS, HCS+, H 2 CS+, 
GH A , CO, NH Z , N 2 , 0-2, H 2 0, Na, H 2 CO, C0 2 , HCN, HNC, 
HNO, H 2 S, C 2 S, HCS, O, OH, H, CH 3 OH, 
CH 3 , NH 2 , HCO, HS, CH 3 0, CH 2 OH 



heat the whole grain, so fn is the fraction of H-atoms rel- 
ative to the ice mantle, plus the refractory core. Thus the 
threshold abundance of H is given by: 

n g (H) = 

where n g (H) is the number of hydrogen atoms in the 
ice mantle, N g is the (average) number of refractory 
atoms/molecules per grain and - n (i) is the sum of all the 
atoms/molecules in the ice mantle. TDulev fc Williams! (120111 ) 
showed that to raise the grain temperature to ~1000K, 12 
kj/mole of stored energy is required, as verified by experi- 
ment. Obviously, much less energy would be needed if the 
required rise in temperature were smaller (e.g. 100 K would 
probably be enough to initiate radical recombination in ice 
mantles). This quantity defines the required concentration 
of reactive species and so, taking the cons ervative value of 12 
kj/mo le, in our standard model we follow iDulev fc Williams! 
(2011) and adopt a value of 0.05 for fn- 

2.2 Phase II 

There are 34 chemical species in the Phase (II) high-density 
expanding gas as listed in Table [3] The (limited) reaction 
network is shown in Table 2] Both the species and the 
reaction network are highly uncertain and ill-determined. 
We speculate that the three-body reactions (for which we 
assume the third body to be H 2 0, the most abundant 
species in the post-sublimation gas-phase) may yield obvi- 
ous organics (such as C 2 HsOH and (CH 2 OH) 2 ). Of course 
other, unidentified, reaction channels/products may also ex- 
ist. Other (small species) products of the reactions are not 
specified in Table [4] Since clear detections of sulfur-bearing 
species in ice mantles do not exist, we do not include any 
complex species/reactions which incorporate sulfur atoms. 
As H 2 is the dominant gas- phase constituent of the Phase 
II gas, we renormalise all integrations to calculate abun- 
dances relative to H 2 0, rather than H nucleons. 

The physical conditions in this extremely high density, 
rapidly expanding, gas are as described in Paper I, so that 
we again consider a sphere of ice, instantaneously sublimated 
into the gas-phase of initial radius ro and with an initial 
density (no) that is comparable to the solid-state number 
density. This gas freely expands, with spherical symmetry, 
into a vacuum at the sound speed v s - Then, as described 



H 2 0, H, H 2 , CO, CH 4 , CH 3 , NH 3 , NH 2 , OH, H 2 CO, HCO, 
C0 2 , H 2 S, HS, CH 3 OH, CH 3 0, CH 2 OH, NH 2 OH, HCOOH, 
C 2 H 6 , CH 3 NH 2 , CH 3 CH 3 0, C 2 H 5 OH, CH 3 CHO, HCONH 2 , 
(CH 2 OH) 2 , CH 2 OHCHO, CH 2 OHNH 2 , CH 3 OCH 3 0, 
CH 3 OH, CH 3 OCH 2 OH 



Table 4. Reaction network in Phase II 



Reaction 
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in Paper I, if ro is assumed to be comparable to the typical 
thickness of an ice mantle and v = 10 4 cms -1 , the evolution 
of the number density is given by 

n/no = 1/ (l + 10 9 i) 3 (1) 

where t is the time in seconds. At such high densities the 
chemistry will be completely dominated by three-body reac- 
tions. Some of the mantle explosion energy will be released 
as heat and we assume that the gas has an initial temper- 
ature of ~ 1000 K, which may help to drive the chemistry. 
We follow the practice of Paper I and adopt a single value 
for the rate coefficient (k 3 s) for all reactions which incorpo- 
rates any implicit dependence on temperature. The radical- 
radical-H 2 reactions are expected to be faster than the 
radical-neutral-H 2 reactions and so we adopt larger values 
for k 3 s than we did in Paper I. As we do not include any 
specific temperature-dependences of the rates, the variation 
of temperature with time is not a relevant parameter in our 
models. 

We follow the chemistry in this expanding sublimated 
gas - which is quite distinct from the background chemistry 
in the surrounding (Phase I) cloud - until quasi-equilibrium 
and/or suppression of chemical activity due to geometrical 
expansion occurs. In practice, this is extremely fast - of the 
order of 1-10 ns (see Paper I); many orders of magnitude 
faster than the chemical timescale of Phase I. This neces- 
sitates a separate numerical integration and the complete 
switching of the chemistry and differential equations at each 
transition between Phases I and II. 
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2.3 Cycling 

After 10 ns, when the Phase II chemical calculations are ter- 
minated, the abundances (calculated relative to H2O) are 
re-normalised to hydrogen and (fully) mixed back into the 
Phase I gas. The cycle, as described above and the duration 
of which is calculated self-consistently, is then repeated. For 
the standard model (£0 = 1.3 x 10 -17 s -1 ) the Phase I dura- 
tion is ~ 3.2 x 10 years and we follow the chemical evolu- 
tion through 5 cycles. We use this number of cycles for two 
reasons: (i) The results indicate that, by then, the solution 
is approaching a limit cycle, and (ii) The total duration of 
~ 1.5 Myr is comparable to the average lifetime for a dense 
cloud. However, for (0 = 1.3 x 10~ 16 s _1 , nn ~ 10 cm" 3 
and so the inter-explosion period is reduced accordingly to 
~ 3.2 — 4 x 10 4 years. We therefore run the calculation for 
50 cycles, so as cover approximately the same time interval 
(~ 1.5 Myr) as for the standard model. 

There will obviously be some optimal value for nn\ if 
it is too high, then that implies that the molecular content 
of the gas is low and ices will not form. If it is too low, 
then the inter-explosion timescale becomes too long for the 
processes described here to be important. The importance of 
tih in controlling mantle growth in molecular clouds may be 
relevant to the determination of the observed critical visual 
extinction for the onset of ice formation. We shall return to 
this question in a later publication. 

Strictly speaking, we should model the gas-phase chem- 
istry and freeze-out of the larger organic molecules created 
in Phase II through successive cycles to ever increasing com- 
plexity. However, we note that the chemistry is very much 
faster and more efficient in Phase II, so we make the sim- 
plifying assumption that the larger species (not included in 
the species list for Phase I) are 'sinked' and simply build up 
in successive cycles, not taking part in the gas-phase chem- 
istry, and not being converted into larger species. By con- 
trast, simple saturated species (such as CH4 and NH3) are 
allowed to re-cycle with a full gas-phase chemistry, freeze- 
out and sublimation. 



3 RESULTS 

At first sight it may seem that there are an unmanageably 
large number of free parameters in the model, but in real- 
ity many of these are degenerate and in effect the chemical 
complexity will largely be defined by two ratios: 

(i) In Phase I; the ratio of the accretion to the inter- 
explosion timescales, and 

(ii) In Phase II; the ratio of the three-body chemistry to 
the expansion timescales 

In the discussion below, we therefore only investigate the 
sensitivities to four parameters; the density in Phase I (nr), 
the cosmic ray ionization rate (Co), the density in Phase II 
(nn) and the rate of formation of radicals in the ice mantle 
(Rrad)- Other parameters, such as the assumed universal 
three-body reaction rate (fcss), have obvious linear effects on 
the large species formation efficiencies (subject to saturation 
limits) . 

The model described above is complex; yielding abun- 
dances as a function of time for both the quiescent phase 



(Phase I) and the sublimation phase (Phase II) for each of 
the cycles of quiescence and explosion. We present the re- 
sults in two forms: (i) as graphical representations of the 
time-dependence of (key) species in the two phases - results 
from all cycles being shown on one plot, and (ii) as tabulated 
abundances of species in Phase I. 

Figures [T] and [2] show the abundances of selected species 
in Phases I and II respectively. The different curves in each 
frame depict the results from different cycles. The time- 
dependences in Phase II, shown in Fig [5] are simple and as 
expected; molecular abundances build up very rapidly until 
saturation (exhaustion of the reactants) and/or geometri- 
cal dilution become important (typically within Ins). From 
Fig [1] it is immediately apparent that, in Phase I, there is 
considerable variation of the abundances with time, partly 
due to the gas-phase chemistry, but mainly as a result of 
the freeze-out of species onto dust grains. Interstellar clouds 
will consist of an incoherent ensemble of material in various 
stages of cycling between the gas-phase and the frozen-out 
states. Bearing these two facts in mind it is therefore appar- 
ent that to make a sensible comparison with observations 
we need to calculate the time- averaged values of the abun- 
dances. For those species where we calculate the full time- 
dependence of the abundances the time-averaged values can 
be easily obtained. For the larger species produced in Phase 
II, whose gas-phase chemistry we do not follow, we use a 
simple analytical approximation. 

On the asumption that freeze-out is the only process 
that affects the abundances of these species in Phase I, the 
rate of change of abundance of species i is simply; 

where ki is the freeze-out rate for species i and n is the den- 
sity. Simple integration then yields, for the average fractional 
abundance; 

(^) = ^#[l-e- fc " lt l 
kiTit L J 

where Xi(0) is the fractional abundance of species i at the 
beginning of Phase I and t is the duration of the phase. 
This approach therefore makes allowance for the statistical 
ensemble of explosion cycles, but does not include the ef- 
fects of sequential enrichment of the gas through successive 
explosions. Whilst some species will approach a limit cycle, 
others will be constantly supplied by successive explosions 
and a cumulative build-up of the explosion products will 
occur. To quantify these effects we have studied the (time- 
averaged) abundances in a sequence of cycles - see Table [(J] 
below. 

We also note that other destruction processes may need 
to be considered. In the case of the larger organic molecules, 
protonation followed by dissociati ve recombination m ay be 
an effective loss mechanism (e.g. iBuckle et al.l [200r3 ). We 
can make a simple, empirically-based, estimate for the ef- 
ficiency of this process as follows: Assume physical param- 
eters as for our standard model and a fractional ionization 
of ~ 10" 8 — 10 -7 , as appropriate for a molecular cloud. As 
an upper limit, if we further assume that the abundance of 
the protonating reactant is equal to the ionization fraction 
and that the reaction rate coefficient is ~ 10~ 9 cm 3 s -1 then 
the upper limit to the protonation (and loss) rate implies a 
destruction timescale of ~ 3 x 10 4 — 3 x 10 5 years. To check 
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this approximation we have run a simple model of a static 
molecular cloud with a full gas-phase chemistry. The ini- 
tial abundance of CH3OH is set high (~ 10 -7 ). The model 
confirms that the CH3OH decays on these timescales. 

These values are comparable to the time interval be- 
tween explosions. By comparison, in our standard model, 
the freeze-out timescale for CH3OH is ~ 2 x 10 s years. 
We therefore conclude that chemical destruction of the 
larger molecules is possibly significant, inhibiting the (time- 
averaged) abundances by a factor of ~ 2 — 5, although this 
effect would be reduced significantly if the period betweeen 
explosions were shorter. 

We present the results from 10 model runs. Model 1 is 
the standard model, which uses the parameter values given 
in Table [1] The variations to these parameter values used in 
the other models are given in Table [S] 

Tables[6][7]and[8]give the (Phase I) time-averaged abun- 
dances for selected species calculated as described above. 
Note that the only species in these tables that are not cre- 
ated exclusively in the explosions are NH3, H2O, H2CO and 
H2S. These are all formed in both the gas and the solid- 
state. In the absence of explosions, all other species would 
have a zero gas-phase abundance and, since we suppress con- 
tinuous desorption in our models, the smaller species would 
also just freeze-out within lMyr. Table [6] shows results for 
the standard model (model 1) in the form of the (average) 
abundances of selected species in each of the five cycles. 
Several conclusions can be drawn from this table: The small 
(Phase I) gas-phase species approach a limit cycle - with 
only small cycle-to-cycle variations after 4 or 5 cycles. The 
same applies to the smaller species produced in the explo- 
sion chemistry (eg. CH3OH, C2H6) which reach abundances 
of ~ 1(T 9 - 1CT 7 . The larger species (eg. CH 2 OHNH 2 , 
C2H5OH) have smaller abundances (~ 1CT 13 - lCT 10 ) but 
these steadily grow as each cycle adds more to the 'reservoir' 
of the larger species. 

Tables[7]and[8]show the sensitivity of the results to vari- 
ations in the free parameters. The results from models 2-4 
are a little surprising at first sight; increasing the (Phase 
I) density drives a faster and more efficient chemistry, so 
that the abundances are somewhat enhanced for model 2 
(n(I) = 10 5 cm -3 ). However at higher densities (in models 
3 and 4), as might be applicable to hot core environments, 
the abundances are notably suppressed - especially for the 
larger species; model 4 yields the lowest abundances of any 
of the models. This is simply a result of the rapid freeze-out 
at high densities so that, although the abundances at the 
beginnning of Phase I are higher, the time-averaged abun- 
dances in Phase I are much lower. Increasing the cosmic ray 
ionization rate (Q - model 5 - again yields relatively small 
enhancements relative to model 1 (by factors of ^ 2 — 6) for 
most species, whilst some are inhibited. However, increasing 
both £ and ni (model 6) results in substantial abundance en- 
hancements, especially for the larger species which in some 
cases (eg. CH3OCH3O) are by more than three orders of 
magnitude. Part of the explanation for this is that higher 
ionization rates result in larger atomic hydrogen abundances 
and hence shorter inter-explosion periods. This, in turn, 
means that the gas-phase depletion factors are lower. 

Models 7 and 8 investigate the effects of a higher den- 
sity in Phase II (nrr). In model 8 the density in Phase I is 
also increased. In the case of model 7 the chemistry of the 



Table 5. Parameter values used in the models. 



Model 


Parameter values 




1 


standard 




2 


n/=10 5 cm -3 , A v =10 




3 


m=10 6 cm -3 , A„=10 




4 


n/=10 7 cm~ 3 , A„=10 




5 


C=1.3xl0 -16 s" 1 , n Hi0 =10cm- 3 , n cyc . 


=50 


6 


C=1.3xl0 -16 s -1 , n/=10 5 cm -3 , n Hfi = 


=10cm~ 3 , 




A„=10, n C!/c . =50 




7 


njj=10 21 cm -3 




8 


nj=10 5 cm -3 , n n =10 21 cm" 3 , A„=10 




9 


R rad =0.02 Myr- 1 




10 


R rad =0.05 Myr" 1 





Table 8. Time-averaged fractional abundances of selected species 
in the final cycle for models 7-10. The nomenclature a(b) implies 
a value of a X 10 6 . 



Species Model 7 Model 8 Model 9 Model 10 



NH 3 


6.3 


-8) 


1 


9( 


-6) 


6 


2 


-8) 


6 


2 


-8) 


H 2 


4.1 


-6) 


1 


6( 


-5) 


i 


1 


-6) 


1 


1 


-6) 


H 2 CO 


1.8 


-9) 


1 


9( 


-8) 


1 


8 


-9) 


1 


8 


-9) 


H 2 S 


6.4 


-10) 


4 


3( 


-9) 


6 


1 


-10) 


6 


3 


-10) 


CH3OH 


2.7 


-7) 


1 


0( 


-7) 


4 


3( 


-7) 


1 


3 


-6) 


NH 2 OH 


5.7 


-8) 


2 


7 


-8) 


8 


4 


-8) 


2 


7 


-7) 


HCOOH 


2.6 


-10) 


1 


3 


-10) 


3 


8 


-10) 


1 


2 


-9) 


C 2 H 6 


5.9 


-8) 


2 


6( 


-8) 


1 


1 


-7) 


2 


9 


-7) 


CH3NH2 


9.; 3 


-9) 


6 


3( 


-9) 


1 


8 


-8) 


4 


7 


-8) 


CH3CH3O 


3.5 


-11) 


5 


4 


-11) 


1 





-10) 


8 


2 


-10) 


C2H5OH 


3.5 


-11) 


5 


4i 


-11) 


1 





-10) 


8 


2 


-10) 


CH3CHO 


1.9 


-11) 


3 


0i 


-11) 


3 


5( 


-11) 


9 


2 


-11) 


HCONH2 


3.8 


-12) 


7 


1( 


-12) 


7 





-12) 


1 


9 


-11) 


(CH 2 OH) 2 


7.0 


-13) 


1 


7( 


-13) 


3 


2 


-12) 


7 


8 


-11) 


CH2OHCHO 


1.5 


-13) 


8 


7( 


-14) 


4 


4 


-13) 


3 


5( 


-12) 


CH2OHNH2 


1.3 


-11) 


1 


3( 


-11) 


3 


9 


-11) 


3 


2 


-10) 


CH3OCH3O 


7.0 


-13) 


1 


7 


-13) 


3 


2 


-12) 


7 


8 


-11) 


CH3OOH 


9.3 


-10) 


1 


9( 


-10) 


2 


1 


-9) 


2 


1 


-8) 


CH3OCH2OH 


6.0 


-13) 


1 


3 


-13) 


2 


7 


-12) 


6 


7 


-11) 



Phase I species is barely affected, whilst the species formed 
in the explosions are enhanced by small factors (typically 
^ 2). Model 8 exhibits a similar behaviour to models 2-4 
in that the higher Phase I density and shorter freeze-out 
timescales results in lower abundances. Models 9 & 10 in- 
vestigate the effects of increasing the radical formation rate 
in the ice mantles (R ra d)- In these cases the abundances of 
both of the radical reactants in the Phase II chemistry are 
enhanced so that there is a second-order effect. As a result 
the abundances of the larger species in model 10 are strongly 
enhanced, by factors of ~ 6 — 400. 



4 DISCUSSION AND CONCLUSIONS 

The results show that, despite the large number of poorly- 
constrained free parameters, the results are fairly insensitive 
within the range of values that we have investigated. Signif- 
icant enhancements of the abundances of large molecular 
species are obtained in the cases of a combination of an el- 
evated cosmic ray ionization together with a high gas-phase 
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Figure 1. Logarithmic fractional abundances (relative to H-nucleons) of selected species as a function of time in Phase I. Results are 
shown, overlaid, for each of the five cycles: 1-solid line, 2-dashcd line, 3-dot-dashed line, 4-wide dotted line, 5-close-dotted line. 



density, or when then the rate of formation of radicals in the 
ices is increased. 

It is interesting to compare the predictions of our dark 
cloud model with observations. However, in so doing, it is 
important to note that a direct comparison is not possible; 
a major conclusion of our study and result from the model 
is that we predict the presence of appreciable abundances 



of large molecules in dark clouds. So far, it has only been 
possible to detect these species towards hot core sources. 
Never the less, despite the very different physical conditions, 
our models are able to reproduce the hot core abundances 
for all species except for the largest ones, e.g. acetamide 
(HCONH 2 ) and glycolaldehyde (CH 2 OHCHO): However we 
note that acetamide has only been observed towards the 
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Figure 2. Logarithmic fractional abundances (relative to H2O) of selected species as a function of time in Phase II. Results are shown, 
overlaid, for each of the five cycles: 1-solid line, 2-dashed line, 3-dot-dashed line, 4- wide-dotted line, 5-close-dotted line. 



Galactic Centre whilst glycolaldehyde has been detected to- 
wards the Galactic Center and one hot core only, hence the 
observed abundances may not be typical. 

Our model is general and not specific to any one source. 
However, although we defer specific modelling and a detailed 
comparison with the dark cloud source TMC-1 (CP) to a fu- 
ture study, we can make some rough comparisons. Some of 



the species given in these tables (NH3, H2O, H2CO, H2S, 
CH3OH, HCOOH and CH3CHO) have been observed in 
T MC-1 (CP) and their abundanc es have been determined 
by ISmith. Herbst fe Chanel |2004) . These abundances are 
generally in very good agreement with our models, with 
the notable exception of methanol (CH3OH) which in most 
of our models is apparently over-produced by a factor of 
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Table 6. Time-averaged fractional abundances of selected species in each cycle for 
the standard model (run 1). A dash indicates a zero or negligibly small abundance. 
The nomenclature a(b) implies a value of a X lO 6 . 



Species Cycle 1 Cycle 2 Cycle 3 Cycle 4 Cycle 5 



NH 3 


1.2(-7) 


2.3(-7) 


5.7( 


-8) 


6.3( 


-8) 


6 


3( 


-8) 


H 2 


1.3(-7) 


1.0(-5) 


4.6( 


-6) 


4.2( 


-6) 


1 


1( 


-6) 


H 2 CO 


4.1(-9) 


7.4(-8) 


7.1( 


-9) 


2.3( 


-9) 


1 


9( 


-9) 


H 2 S 


9.5(-14) 


7.4(-10) 


6.3( 


-10) 


6.4( 


-10) 


6 


4( 


-10) 


CH3OH 




1.4(-7) 


1.5( 


i-t\ 

-') 


1.6( 


-•I 


1 


6( 


-') 


NH 2 OH 




2.0(-8) 


2.4( 


-8) 


2.6( 


-8) 


2 


9( 


-8) 


HCOOH 




8.2(-12) 


1.1( 


-10) 


1.2( 


-10) 


1 


3( 


-10) 


C 2 H 6 




4.7(-8) 


4.9( 


-8) 


4.9( 


-8) 


4 


9( 


-8) 


CH3NH2 




7.0(-9) 


7.3( 


-9) 


7.4( 


-9) 


7 


5( 


-9) 


CH3CH3O 






8.8( 


-12) 


1.2( 


-11) 


1 


5( 


-11) 


C 2 H 5 OH 






8.8( 


-12) 


1.2( 


-11) 


1 


5( 


-11) 


CH3CHO 




2.8(-12) 


1.4( 


-11) 


1.4( 


-11) 


1 


1( 


-11) 


HCONH 2 




4.2(-13) 


2.5( 


-12) 


2.7( 


-12) 


2 


8( 


-12) 


(CH 2 OH) 2 






4.5( 


-14) 


1.1( 


-13) 


1 


8( 


-13) 


CH 2 OHCHO 






5.6( 


-14) 


6.3( 


-14) 


6 


5( 


-14) 


CH 2 OHNH 2 






1.7( 


-12) 


3.7( 


-12) 


5 


7( 


-12) 


CH 3 OCH 3 






4.5( 


-14) 


1.1( 


-13) 


1 


8( 


-13) 


CH3OOH 






7.8( 


-11) 


1.7( 


-10) 


2 


6( 


-10) 


CH 3 OCH 2 OH 






3.9( 


-14) 


9.3( 


-14) 


1 


5( 


-13) 



Table 7. Time-averaged fractional abundances of selected species in the final cycle for models 
1-6. The nomenclature a(b) implies a value of a X 10 . 



Species Model 1 Model 2 Model 3 Model 4 Model 5 Model 6 



NH 3 


6.3( 


-8) 


1.9( 


-6) 


3.3 


-7) 


3.5 


-8) 


1.1( 


-7) 


5.0( 


-8) 


H 2 


4.1( 


-6) 


1.6( 


-5) 


1.6 


-6) 


1.6 


-7) 


7.3 


-6) 


1.2 


-5) 


H 2 CO 


1.9 


-9) 


2.5( 


-8) 


4.5 


-10) 


2.4 


-11) 


1.0 


-9) 


1.8( 


-9) 


H 2 S 


6.4( 


-10) 


4.3( 


-9) 


5.2 


-10) 


5.2 


-11) 


7.9( 


-10) 


4.3( 


-9) 


CH 3 OH 


1.6( 


-7) 


5.5( 


-8) 


5.5 


-9) 


5.5 


-10) 


3.2 


-7) 


8.7( 


-7) 


NH 2 OH 


2.9( 


-8) 


1.3 


-8) 


1.7 


-9) 


1.7 


-10) 


9.2( 


-8) 


3.9 


-8) 


HCOOH 


1.3 


-10) 


7.7 


-11) 


1.3 


-12) 


6.9 


-14) 


4.7( 


-11) 


4.7( 


-9) 


C 2 Hg 


4.9( 


-8) 


2.1( 


-8) 


2.1 


-9) 


2.1 


-10) 


9.9( 


-8) 


2.2( 


-7) 


CH 3 NH 2 


7.5( 


-9) 


4.7( 


-9) 


6.3 


-10) 


6.5 


-11) 


1.3 


-8) 


1.2 


-8) 


CH3CH3O 


1.5 


-11) 


2.2 


-11) 


2.3 


-12) 


2.2 


-13) 


4.4( 


-11) 


2.0( 


-9) 


C 2 H 5 OH 


1.5 


-11) 


2.2( 


-11) 


2.3 


-12) 


2.2 


-13) 


4.4( 


-11) 


2.0( 


-9) 


CH3CHO 


1.4 


-11) 


2.9( 


-11) 


4.8 


-13) 


2.6 


-14) 


6.6 


-12) 


1.1 


-9) 


HCONH 2 


2.8( 


-12) 


6.3 


-12) 


1.4 


-13) 


8.1 


-15) 


9.1( 


-13) 


1.9 


-11) 


(CH 2 OH) 2 


1.8 


-13) 


3.6( 


-14) 


3.5 


-15) 


3.5 


-16) 


8.5( 


-13) 


2.1( 


-10) 


CH 2 OHCHO 


6.5( 


-14) 


4.3( 


-14) 


7.1 


-16) 


3.7 


-17) 


1.6 


-14) 


6.8( 


-12) 


CH 2 OHNH 2 


5.7( 


-12) 


4.9 


-12) 


6.7 


-13) 


6.9 


-14) 


3.3 


-11) 


7.2( 


-11) 


CH3OCH3O 


1.8( 


-13) 


3.6( 


-14) 


3.5 


-15) 


3.5 


-16) 


8.5( 


-13) 


2.1( 


-10) 


CH3OOH 


2.6( 


-10) 


5.1 


-11) 


5.1 


-12) 


5.0 


-13) 


1.6 


-9) 


2.1( 


-8) 


CH 3 OCH 2 OH 


1.5 


-13) 


2.7( 


-14) 


2.7 


-15) 


2.6 


-16) 


8.3( 


-13) 


1.8 


-10) 



up to ~100. However, there are two reasons why this is 
probably not discrepant; (i) there is considerable flexibil- 
ity/uncertainty in the physical parameters (given in Ta- 
ble [TJ which determine the efficiency of the formation of 
complex molecules, and (ii) as described in the previous sec- 
tion, we have probably underestimated the loss rate for large 
molecules in the quiescent phase (Phase I). Our model also 
over-produces water. Howe ver, the clumps i n TMC-1 are 
known to have substructure l|Peng et al.lll99cj ) which means 
that ultraviolet penetration into CP is easier than it would 
be in a uniform region - as we have assumed. So, we would 
expect that the time-averaged H 2 abundance should be 



lower than the values we have computed for the non-specific 
case. In addition, we also note that the efficiency of the 
surface reaction converting O to H 2 is required by obser- 
vations to be high, but is undetermined. We have assumed 
a value of 100% but a factor of a few less than this would 
still be compatible with the observed strengths of the water 
ice features. Looking ahead, we expect that the predicted 
presence of hot-core type molecules within molecular clouds 
are at levels that can be detected with new facilities, such 
as ALMA. 

As it stands, the model is largely hypothetical. We have 
been able to give reasonable estimates for the key parame- 
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ters. However, we require further experimental study of the 
explosion mechanism. In addition, we need confirmation of 
the viability of the proposed three-body reaction pathways 
and quantification of the reaction rates. However, several 
conclusions can already be drawn from this study; 

• From our results it can be seen that our model pre- 
dicts the presence of detectable (observationally verifiable) 
gas-phase abundances of 'large species' in molecular clouds. 
However, the mechanism that we describe has a limiting ef- 
fect on the size of molecules that can be made, because of 
the very short time available for rapid 3-body reactions to 
take place. 

• Our model also suggest that there could be significant 
amounts of undetectable molecules, such as C2H6, present 
in dark clouds. If so, these species - which tend to be ignored 
in chemical networks - may also play a role in the formation 
of observed larger species, through ordinary gas phase ion- 
molecule reactions. For example, a reaction sequence involv- 
ing C2H6 and CH/C2H radicals can lead to t he formation 
of be nzene (CeHg) in dark cloud environments (I Jones et al.l 

• In recent years, spectral line surveys of dark cores, such 
as those of TMC-1, have revealed a surprising repertoire 
of complex molecules and large carbon chain species includ- 
ing propylene (CH2CHCH3). This latter species, expected in 
hot core s rather than in dark clouds, has been detected in 
TMC-1 (|Marcelino et al.ll2007r i but not in Orion. In general, 
other more common species, such as HNCO and CH3OH, 
are suprisingly abundant in some dark cores. Pure gas-phase 
models at low temperatures (~10 K) cannot reproduce their 
abundances and grain surface reactions on icy mantles are 
often invoked. However, it is not clear what mechanisms are 
responsible for the release of icy mantles in cold cores; it is 
unlikely that non-thermal desorption alone would suffice and 
the presence of shocks inducing explosive inj ection of grai n 
mantles has been suggested in the past (e.g. ISnvder]l2006|) . 
Our model provides a novel mechanism for the production 
of large complex molecules in cold dark clouds. 

• Methanol (CH3OH) and formaldehyde (H 2 CO) are 
found in ices in quiescent regions. In our model these species 
are synthesised in the gas phase following the ice mantle ex- 
plosions, are ejected and subsequently frozen onto the grains 
in the next cycle. This is different from the conventional pic- 
ture in which they are made by the hydrogenation of CO 
in the ice. Although we allow for the partial conversion of 
CO to CO2 we have deliberately suppressed the solid-state 
conversion of CO to CH3OH in our models, so as to clearly 
separate the significance of our mechanism from the conven- 
tional channel. We find that both processes may be effective. 

• The composition of the ice mantles and the nature 
of gas-grain interactions differ from previous models. We 
have postulated an alternative desorption mechanism which 
should be considered alongside other desorption processes. 
Our model predicts that the ices will probably be more 
chemically complex than previously recognised. Moreover, 
the chemical processing is cumulative, so that we expect the 
the chemical complexity of ices in molecular clouds to grow 
in time. 

• The presence of trace metal elements, such as Na, in the 
hot sublimated gas together with OH radicals may lead to 
the formation of NaOH. NaOH has a large dipole moment, a 



rotational constant of Bo ~ 12.57 GH z and an observation- 
ally detectable rotational spectrum l|Pearson fc Truebloodl 

• The chemical network that we employ in the high den- 
sity phase (Phase II) is highly speculative, probably incom- 
plete, and somewhat pessimistic in that it only involves re- 
actions with H2O as a passive third body; i.e. one that stabi- 
lizes associating species, but does not chemically react with 
them. Our results should therefore be considered as lower 
limits to the large molecule formation efficiencies. We can 
speculate, for example, that the abundant species CO (which 
is technically a radical) could be involved as an active third 
body in which it chemically combines with the other reac- 
tants. If, for example, these were NH2 and CH2OH radicals 
then it may be possible to form a large molecule, such as 
glycine (NH2CH2COOH), in a single step process. 

• The model may provide a mechanism for the for- 
mation of larger molecules, of biochemical importance, in 
molecular clouds. If we consider C2H4O2, then three iso- 
mers of this molecule are detected in hot cores; acetic 
acid (CH3COOH), methyl formate (HCOOCH3) and gly- 
colaldehyde (CH2OHCHO). No detections of these species 
have been made in molecular clouds. However, in hot 
cores, acetic acid and glycolaldehyde are only detected 
rarely (such as in Sgr B2(N) and W51), whilst methyl for- 
mate is more abundant and is observed in many sources. 
Thus, for example, the relative abundances of (acetic 
acid): (glycolaldehyde): (methyl formate) are ~1:4:26 in Sgr 
B2(N). This is consistent with the observation that interstel- 
lar molecules with a C-O-C backbone structure are preferred 
over those with a C-C-O structure (|Hollis. Lovas fc Jewelll 
l200d) . 

Our model only specifically identifies glycolaldehyde as a 
reaction product. It is possible that methyl formate could be 
formed as a result of a three-body reaction involving the rad- 
icals CH3O and HCO. However, even if this assumption is 
made, then the model predicts methyl formate abundances 
that are less than those of glycolaldehyde. There are three 
possible explanations for this apparent discrepancy: (i) The 
physical and chemical conditions in hot cores are very dif- 
ferent to molecular clouds, so that molecular ratios deter- 
mined from hot core observations may not apply, (ii) as is 
highly likely to be the case, our Phase II reaction network is 
neither fully accurate, nor complete and/or (iii) alternative 
mechanisms may be operating which produce the methyl for- 
mate. Indeed a variety of mechanisms have been proposed 
for methyl formate formation, including the solid-state re- 
action of methanol and CO in cosmic-ray irradiated ices 
(|Modica fc Palumbo I l20ld ; lOcchiogrosso et all l201ll ). the 
solid-state reaction of formic acid and methanol (Hollis et al. 
2000) and ga s-phase reaction of methanol and f ormaldehyde 
in hot cores (|Caselli. Hasegawa fc HerbstHl993D . 

Looking at the structure of these three isomers we see that 
methyl formate, glycolaldehde and acetic acid have a cen- 
tral O-atom (ether group), methylene group and carbonyl 
group respectively. As we do not have atomic oxygen in our 
Phase II chemistry it is therefore perhaps not surprising that 
methyl formate is not produced in our model. 

If we go on to speculate about the possible formation of 
larger bio-molecules then we should note that methyl for- 
mate is less important that the other two (C-C-O structure) 
isomers; glycolaldehyde is the simplest sugar, whilst acetic 
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acid is only an NH2 group away from the simplest amino 
acid, glycine (NH 2 CH 2 COOH). Indeed, it is highly likely 
that the main synthesis channels for glycine involve acetic 
acid so that it may be so mething of an essential pre-cursor 
to bio- molecule formation l|Snvderll2006I i. 

Whilst acetic acid is not specifically included in our reac- 
tion scheme, it is a perfectly feasible product of the reactants 
that are present in the high density phase (Phase II) of our 
model - indeed more so than methyl formate. It is also worth 
noting that our model does produce formic acid (HCOOH) 
with typical fractional abundances of ~ 10 -10 — 5 x 10 -9 . 
Formic acid is al so believed to be a significant species in 
glycine synthesis |Snvderl l 2006). An alternative mechanism 
for the formation of complex organic molecules was proposed 
bv lSorrelll i|200ll ) in which radicals are created in ice mantles 
by UV photolysis. In that model, the grains are heated by 
grain-grain collisions and the radicals react with each other 
in the mantles. Chemical explosions release the complex or- 
ganics, including glycine and glycolaldehyde into the gas- 
phase. There is no reas on why the mechanism that we are 
propo sing and that of iGarrod. Widicus Weaver fc HerbstI 
(2008) should not be both possible and may indeed com- 
plement each other. 
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